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ABSTRACT 

Five observations of the hard X-ray spectrum of Her X-1 from RXTE show 
that the ~ 41 keV energy of the cyclotron scattering line is constant within 
statistics of a few percent per observation. The overall spectral shape, on 
the other hand, varies somewhat, with an RMS of 2%. If the 41 keV feature 
truly originates as cyclotron resonance scattering in an unchanging 3 x 10^^ 
Gauss dipole held not far above the neutron star surface, these observations 
constrain the average height of scattering to within a range of 180 meters. This 
is consistent with models which put the radiating structure within meters of 
the surface of the neutron star. In other pulsars observed line centroid changes 
have been correlated with luminosity changes, and if interpreted as variations 
of the height at which scattering takes place, many hundreds of meters are 
required. These RXTE data, which sample nearly a factor of two in unabsorbed 
luminosity, are in conflict with a particular model for such an extended radiating 
structure. Comparison with other observations over many years indicates 
strongly that the centroid energy of this absorption line has increased some time 
between 1991 and 1993 by 23%, from 34 keV to 41 keV. Moreover, the cutoff 
energy of the spectral continuum increased at the same time from 16 keV to 
20 keV, which is, within the statistical error of 5%, in direct proportion to the 
centroid. This may be a sign that both these characteristics of the spectrum 
are controlled in the same way by the magnetic held strength in the region of 
scattering. 


Subject headings: pulsars: individual (Her X-1) — X-rays:stars 



1. Introduction 


The accreting X-ray pulsar, Hercules X-1, emits a spectrum with a localized depression 
at roughly 40 keV which is generally interpreted as an absorption line resulting from 
scattering of photons on electrons whose motions are constrained in one spatial dimension 
to transitions between Landau levels in a teraGauss magnetic held at the neutron star’s 
polar cap. Commonly referred to as a cyclotron line, this feature was discovered in 1977 
(Triimper et ah 1978), although hrst identihed phenomenologically as an emission feature 
at somewhat higher energy. This feature in Her X-1 has since been observed repeatedly 
(e.g. Gruber et al. 1980, Voges 1984, Tueller et al. 1984, Soong et al. 1988, Mihara 1995, 
Kunz 1995, Dal Fiume et al. 1998). Similar features have also been identihed in the hard 
X-ray spectra of many of the accreting X-ray pulsars (e.g. Mihara 1995, Heindl et al. 1999, 
Dal Fiume et al. 2000). Cyclotron lines ahord a diagnostic of plasmas under extreme 
conditions and can provide a key to the structure of the magnetosphere and the mass how. 

Mihara et al. (1998) have reported a study of repeated observations with Ginga of 
cyclotron lines in a number of X-ray pulsars. In a few instances, although not with Her X-1, 
they found a correlation with X-ray luminosity which was consistent with the cyclotron 
scattering taking place from hundreds to many hundreds of meters above the neutron star 
surface, at the peak of an “accretion mound” structure proposed by Burnard, Arons & Klein 
(BAK, 1991). BAK predict a strong dependence of this height on accretion how. Thus in 
this model variability of the cyclotron line energy can also be an important diagnostic of 
the scattering region and accretion process. 

Studies of line formation (Meszaros & Nagel 1985, Isenberg, Lamb & Wang 1998, 
Araya-Gochez & Harding 2000) have indicated the possibility of line prohles which deviate 
from the simple Gaussian and Lorentzian prohles which have been used for analysis to 
date. Such deviations are evident as particularly broad lines (Heindl 1999, Orlandini 1999). 
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However, the line of Her X-1 is narrower, and measnrement of strnctnre beyond the hrst 
and second moments (centroid and width) reqnires detector resolntion higher than the 
5E/E of a few obtainable with inorganic scintillators, as well as very large collecting area 
and observing time. 

Using data from hve observations over two years with the large-area instrnments 
aboard Rossi X-Ray Timing Explorer (RXTE) we have performed detailed spectroscopy 
of the spectral continnnm, the cyclotron line and the flnorescent iron line of Her X-1. We 
hnd only minor variability between these observations. Exercising carefnl data selection, 
we have compared the resnlts with those from earlier measnrements between 1977 and 
1994. Eor this pnrpose the HEAO-1 data were freshly analyzed with modern software and 
spectral models. To a level of perhaps 10% the pre-1991 spectral shapes are consistent with 
each other, as are the spectra from 1994 onwards. Between 1991 and 1994 there strongly 
appears to have been a change of the cyclotron line centroid and the cntoff energy of the 
continnnm. This is in marked contrast to the power law index and exponential fold energy, 
for which the average valnes remained constant to within a few percent. In §2 we describe 
instrnmentation and observations, in §3 analysis and resnlts, and in §4 we briefly discuss 
possible interpretations. Preliminary results from a portion of this RXTE data set have 
been reported earlier (Gruber et ah 1997, 1998, 1999). 


2. Instruments and Observations 

The RXTE, launched late in 1995 into low earth orbit, performs pointed observations 
with the Proportional Counter Array (PGA, Jahoda et ah 1996) and the High Energy 
X-Ray Timing Experiment (HEXTE, Rothschild et al. 1998). The PGA consists of hve 
Xenon proportional counters (called PCU’s) sensitive in the energy range 2-60 keV with 
1400 cm^ each. The HEXTE consists of eight Nal scintillator detectors of area 200 cm^ each 
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and energy range 15-250 keV. Both instruments are non-imaging and have coaligned fields 
of view collimated to 1° FWHM. Time resolution was microseconds. Energy resolution was 
16% at 6 keV for PCA, and 15% at 60 keV for HEXTE. 

The five RXTE observations (Table 1) each consisted of several satellite orbits, each 
orbit typically containing 3300 seconds on-source. All five PCU’s were used except for the 
observations of 1996 August, for which only four detectors were operated. After the first 
observation one of the eight HEXTE detectors was lost for spectroscopic purposes. The 
observations were collected over two years from 1996 and 1997. Each was selected to be 
near the peak of the Her X-1 “main-on” state and to avoid periods of eclipse and other 
forms of obscuration, such as pre-eclipse dips. Exposures were varied, but averaged 15000 s 
per observation. 

The study of long-term changes relies particularly on archival HEAO-1 A-4 (Matteson 
1978) data, which we re-analyzed with modern software and spectral models. The two 
Low-Energy Detectors of this experiment each had area 103 cm^ and spectral resolution 
25% at 60 keV. The A-4 instrument was quite accurately calibrated. The calibration is 
described in part in Jung (1986) and Gruber et ah (1989), and additional relevant details 
are given in the Appendix. 


3. Analysis and Results 

The spectral form characteristic (White, Swank & Holt 1982) of accreting X-ray 
pulsars, a power law of photon number index P, E“^, and with high-energy exponential 
rolloff Ef above a cutoff energy Ec, , was the starting point for analysis. To 

represent cyclotron absorption we employed the multiplicative profile 

^(-T-exp{-{E-Erf /2aA) ^ 
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where is the centroid, a is the line RMS, and r is the optical depth. Earlier work in 
the literature began with additive and subtractive forms to represent line emission. With 
the general acceptance of the line formation process as resonance cyclotron scattering 
of photons from the continuum, multiplicative prohles were adopted. With Ginga and 
BeppoSAX a nearly Lorentzian prohle was also employed: 

^{-T-{a-EIEr)l{(E-Erf+a^)^ 

with r, Er and a as above. We report detailed results only with the Gaussian prohle in 
optical depth, but note that best-ht values for the line centroid differ only by -5% using 
the Lorentzian prohle, and a similar amount for subtractive prohles. All results of this 
analysis are signihcant also if the Lorentzian prohle is used in place of the Gaussian. We 
note that the employment of any of these prohles is largely historical and heuristic. Recent 
calculations (Araya-Gochez & Harding 2000) indicate a wide variety of possible line prohles. 

Spectra averaged over pulse phase were employed for this work. A major goal was 
comparison with historical observations, many of which are available only in this form. 
Although spectra can be strongly phase dependent in some pulsars, with Her X-1 only 
modest changes with phase have been observed for selected parameters, in particular the 
power law index change of 30% (Pravdo et al. 1977) and a 20% change of the cyclotron 
line centroid (Gruber et al. 1980, Voges 1984). Pulse phase-resolved spectroscopy will be 
reported in a later publication. The observations were collected at various orbital phases 
(Table 1), but we expect no ehect from this, because a lack of spectral changes with orbit 
was reported by Pravdo (1976). Moreover, the eight HEAO-1 observations of Soong et 
al. (1990) showed no orbital variation. We performed spectral analysis with the cut-off 
power law continuum model already mentioned, and we tried and rejected another model 
in current use, the Fermi-Dirac (Tanaka 1986). These forms are used empirically only: we 
are unaware of any theoretical forms for the spectral continuum which permit meaningful 



comparison with the observed spectra. 


3.1. RXTE Spectral Analysis 

Data were selected from the five observations such that the source was uneclipsed 
either by the primary or the earth and such that characteristic absorption events (“dips”, 
Nh > lO^^cm”^) were absent. A substantial dead-time correction was applied to the 
HEXTE data. The HEXTE background was directly measured by switching the HEXTE 
look direction on- and off-source every 16 seconds, while the PCA internal background was 
estimated with the RXTE software tool PCABACKEST and the strong-source background 
model. This model is known to have an accuracy of the order of a few percent. Following 
established procedure (e.g. Wilms et ah 1999) we tuned the background estimate to the 
observed high-energy portion of the data by adjusting a scaling factor to produce zero PCA 
net counts above 60 keV. The RMS correction for the hve observations was 4.4%. As a 
further measure to limit sensitivity to residual background subtraction errors, PCA data 
below 3.5 keV and above 25 keV were not used, except as a check. Within the 3.5 - 25 keV 
range the average source spectrum was at least four times background, and usually much 
greater. Even with this favorable signal-to-background residuals to fits initially showed 
structure at the 1% level. We found very similar structure in power-law hts to the Crab 
spectrum with PCA, and we attributed this to imperfectly modeled background at the level 
of accuracy of the background model, about one milliCrab. We were able to model the 
spectrum of the extra background in the Crab observations as an extremely broad, a = 
2.5 keV, Gaussian line centered at about 5 keV. We then applied it successfully to remove 
the dominant residuals in the hts to Her X-1, and the adjusted constant of normalization 
again yielded a hux of close to one milliCrab, as with the Crab itself. Except for this extra 
background component and for details, our procedure for dealing with these systematic 



errors closely follows that described by Wilms et al. (1999, see Appendix B). As a practical 
matter, the HEXTE data above 25 keV are statistically much stronger than the PCA, 
and very little statistical weight is gained by adding PCA data above 25 keV. As a check, 
however, we have performed joint hts including PCA data to 60 keV. The higher-energy 
PCA data are quite consistent with the HEXTE data. In particular, in each observation 
the PCA independently confirms the presence of the absorption line at about 40 keV. The 
HEXTE data range extended from 20 keV, chosen to be above the electronic threshold for 
all observations, to 100 keV, above which energy data were discarded because the Her X-1 
signal was below detectability. 

Standard PCA energy response matrices produced with version 5.0 ftools were 
employed. These were checked against contemporary Crab observations, and non-statistical 
errors of 0.25% RMS were thus measured. This error was applied to each PCA energy 
channel as a systematic error. This 0.25% error represents a rehnement of the calibration 
procedures which earlier led Heindl et al. (1999) to assign a 1% systematic error. The 
residuals to the Her X-1 hts with this high-energy cutoff model for the continuum showed a 
depression at the cut energy. This was identihed as an artifact of this spectral model, which 
has discontinuous slope at the cut energy. This artifact was quite effectively suppressed by 
modeling a Gaussian absorption at the cut energy with optical depth 0.15 and sigma 2.5 
keV. We note that Burderi et al. (2000), in analyzing BeppoSAX data for Cen X-3, also 
found the need to introduce an empirical rounding function at the cutoff energy with this 
continuum. They used a polynomial rather than a Gaussian form. 

With this model and the very low 0.25% systematic error for PCA, hts to all hve 
observations produced acceptable chi-squares and a hat and featureless pattern of residuals. 
In particular, we found no evidence for the 18 keV spectral break reported from a 
BeppoSAXohservaiion (Dal Fiume et al. 1998). 



9 


Fits to a power law with the Fermi-Dirac (Tanaka 1986) cutoff function, 

also with fold and cutoff energies but with continuous derivative, were tried in the hope 
of modeling the spectral turnover more simply. Instead, a pattern of correlated residuals 
over the entire PCA energy range resulted, strongly indicating a poor £t and functional 
mismatch with the data. This test with an alternative continuum did show, however, that 
best-£t values for the line centroid were even more insensitive to the choice of continuum 
than their 5% sensitivity to choice of line prohle. 

In Figure 1 we show the best £t and residuals for observation 2, which had the longest 
duration. The middle panel shows the large residuals, particularly at 40 keV, of a £t 
without cyclotron line, while the lowest panel shows the acceptable fit when the line is 
included. The displayed PCA data include the 0.25% systematic error per datum. 

In Table 2 we list best-£t parameters and errors from the hve observations for the 
successful continuum model with power law and high-energy cutoff, plus iron K emission 
and cyclotron scattering feature. In the last column of Table 2 are shown results of a 
simultaneous fit to all five observations with free normalizations for each but with the same 
spectral shape parameters, including line depth. For the £t to the combined data the large 
Xr of ~4.6 is dominated by small but real spectral changes from observation to observation, 
primarily in the spectral index and cutoff energy. We estimate the spectrum-to-spectrum 
variability to be 1.8% RMS. 

Best-fit parameter values in Table 2 were determined by minimizing the statistic, 
using the systematic plus statistical errors for PCA spectral data and statistical errors only 
for HEXTE. The parameter error estimates were determined from a number of Monte-Carlo 
realizations using xspec: starting with a best-£t model we determined the average and RMS 
scatter for each parameter. For the parameters which dehne the continuum shape these 



10 


experimentally-determined errors averaged a factor 1.8 greater than the values reported 
directly by xspec, while corresponding estimates for the parameters of the scattering line 
were larger only by a factor 1.1. These estimates are smaller than those obtained by the 
prescription of Lampton, Margon & Bowyer (LMB, 1976), which gave factors averaging 2.2 
larger than the direct xspec estimates. 

Given the satisfactory results obtained with the best-fit spectral model, we could also 
search for variability of the individual parameters of this model. For the hve measured 
values of each parameter in Table 2 we calculated averages and chi-squares. With either 
our Monte-Carlo estimates or the LMB errors large chi-squares for the parameters for the 
continuum shape strongly indicated very signihcant but nevertheless modest variability. 
Conversely, for the scattering line parameters both methods indicated stability, with 
reduced chi-squares near unity or below. Reduced chi-squares using the LMB estimates for 
the three line parameter errors were all about 0.25, indicating (P=0.001) that the LMB 
errors are overestimated. 

The parameters which depend most strongly on the correction for PCA systematic 
error are the power law index and the iron line parameters. The cutoff, fold and cyclotron 
line parameters are determined largely by the HEXTE data, for which pure statistical 
errors were employed. 

Of special interest is the possible variability of the cyclotron line centroid. The hve 
RXTE centroid values are consistent with no variability. The chi-square of 5.7 (4 dof) for 
the centroid values permits an extra random variability as high as 1.5 keV (90% conhdence) 
in the centroid energy. 

The hve RXTE observations were each made in the hrst few days of a Her X-1 35-day 
main-on state. Strictly speaking, then, the measured 1.8% spectral stability applies only to 
the early part of the main-on. The spectrum varies somewhat with pulse phase (Soong et 
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al. 1990) and the pulsation shape is known to vary regularly with 35-day phase (Scott et 
ah 1997). So there is reason to expect spectral variability with 35-day phase. However, 
the HEAO-1 observations (Soong et al. 1990), which sampled the main-on state fairly well, 
showed no change of pulse phase-average spectral shape with 35-day phase. 

The contemporary observation with BeppoSAX (Dal Fiume et al. 1998), which has 
similar statistical weight to the RXTE observations, yielded similar results. The BeppoSAX 
values for the cyclotron line parameters agree with the RXTE values to within the errors 
of Table 2. There is reasonable agreement for the iron line parameters. The continuum 
parameters are harder to compare because of the 18 keV break in the BeppoSAX power 
law index, which we do not conhrm. In the next section we compare these results with 
observations made well before the launch of RXTE. This comparison indicates signihcant 
changes for the cyclotron line energy and the cutoff energy of the continuum. Significant 
change was not found in other parameters of the spectrum. 

The 41 keV absorption line was observed at high significance in all five RXTE 
observations. Triimper et al. (1978), interpreting the cyclotron feature as emission at 
55 keV, reported a possible harmonic overtone, although this has never been conhrmed. 
Multiple harmonic lines have been observed in the spectrum of a few other x-ray pulsars, 
most notably 4U0115-I-63 (Heindl et al. 1999; Santangelo et al. 1999), which has four or 
five. We have searched for a harmonic overtone of the 41 keV scattering feature by means of 
a simultaneous £t to the HEXTE data of all five observations. The resulting best-fit optical 
depth for the overtone at 82 keV is zero, and the 95% confidence upper limit is 0.23. This 
limit shows that the cyclotron scattering is quite different in Her X-1 than in 4U0115-I-63, 
whose overtones have optical depths greater than unity. 



- 12 - 

3.2. Reanalysis of HE A 0-1 

To permit more reliable comparison with earlier work we have re-analyzed archival 
HEAO-1 A2 and A4 spectra using the same high-energy cutoff and line prohle functions 
employed here with RXTE data. Another important reason for re-analysis is that 
considerable work (see Appendix) has gone into the calibration of the HEAO-1 A4 since 
the results of Soong et ah (1990). For two of the HEAO-1 observations we were able to 
£t simultaneously 1-18 keV data from an Argon counter of the A2 experiment (Rothschild 
et ah 1978), thereby obtaining a good measure of the power law index. One such fit is 
shown in Figure 2. Two sets of residuals are shown in Figure 2: the middle panel contains 
residuals to a fit without the cyclotron line and clearly shows a pronounced deficit at 35 
keV; the bottom panel contains acceptable residuals to a fit which includes the line. 

The HEAO-1 spectral parameter values and errors shown in Table 3 are from this 
re-analysis. By comparison with Soong et al. (1990), values differ somewhat. The spectral 
index and cutoff energy are smaller by 20%, the fold energy is 10% larger, the line centroid 
(3% smaller) is nearly the same, and the line width is smaller by 30%. The change of line 
width is probably not significant, given the resolution and statistics. As with the RXTE 
data, we tested the sensitivity of the cyclotron line centroid value to different line profiles, 
as well as the two continuum models. We found changes only of the order of 1 keV, even 
less than with RXTE, Errors were estimated as described here for RXTE. They are a factor 
of 2-3 larger than those reported by Soong et al. (1990), who used a simpler procedure to 
estimate errors. 

The most interesting result from the reanalysis are the HEAO-1 values for the cutoff, 
which lie between 16 and 17 keV, considerably smaller than the 24 keV reported as a 
best-fit for HEAO-1 on-pulse minus off-pulse spectra by Gruber et al. (1980), and also less 
than the 21 keV assumed by Soong et al. (1990). These two authors did not have the use of 
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the HEAO-1 A2 data, which powerfully constrains the value of the spectral index. Indeed, 
the values of the spectral index, cutoff and fold are very highly correlated in the A4 data, 
because these data have a lower threshold at 13 keV. The correlation has the sense that 
overestimating the cntoff energy reqnires the £t to compensate by moving the index to a 
larger value and the fold energy to a smaller valne. When we fit to the A4 data alone, the 
index moved to a valne above nnity, the cntoff energy moved to 18 keV, and the fold energy 
to 13 keV. A cutoff energy at 21 keV, as used by Soong et ah (1990), requires a chi-square 
higher than this A4-only best-fit by ~4. Such a difference is reasonably consistent with 
small changes in the response matrix, backgronnd subtraction, and fnnctional forms from 
those used earlier. We thus believe that the earlier cutoff values of Gruber et ah (1980) 
and Soong et ah (1990) were considerably overestimated, largely dne to the limited spectral 
range of the A4 data. We note that the HEAO-1 A2/A4 data sets are very well matched to 
the RXTE PCA/HEXTE data sets in spectral range and statistical significance. Therefore 
comparisons of spectral parameters between observations performed with these missions is 
much more reliable than between other pairs of missions. 


3.3. Long Term Changes 

The valnes for the line centroid in the present work are ronghly consistent with 
measnrements of valnes near 44 keV obtained between 1993 and 1995 with CGRO/BATSE 
(Freeman et. al 1996). The RXTE, BeppoSAX and BATSE valnes differ from all earlier 
measnrements, which centered near 35 keV (Kunz 1995). Earlier data with good statistics 
and reliable energy calibrations inclnde those from the HEAO-1 A4 experiment (Soong 
1988), the GRIS germanium detector (Tueller et al. 1984), the HEXE experiment aboard 
the Russian platform MIR (Kunz 1995), and Ginga (Mihara 1995). However, functional 
forms for the representation of the continnnm and line profile have been various, making 
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comparison somewhat uncertain. In Table 3 we have collected reported results for continuum 
(power law with a high energy exponential cut off) and cyclotron line parameters. The iron 
line parameters, which did not change sensibly, are ignored. 

The cyclotron energy between 1976 and 1996 is charted in Figure 3. Clearly, the data 
before 1991 are consistent with an average value near 35 keV and the data after 1991 are 
consistent with a higher average of 41 keV. When we analyze the parameters in Table 3 
for the formal signihcance of a change of average value (a step function in time) using 
the F-test, we obtain the average values and significances shown in Table 4. The large 
reduced chi-squares are perhaps inevitable given calibration errors and somewhat differing 
functional forms for hts, as well as the possibility of variability on any time scale. Formal 
statistical results based on such high chi-squares are not valid. However the F-statistic, 
which is based here on the comparison of chi-squares for the null hypothesis of no change, 

and for a change between 1991 and 1994, x}, should be more robust and provide at 
least an upper limit to the the null hypothesis (no change) probability. With this statistic 
the line centroid has in fact a very low signihcance < 1.8 x 10“^^ for no change. Probably 
signihcant is a change of the cutoff energy, with Pnuii < 0.01. Possibly signihcant at 90% 
conhdence are small changes of the power law index and fold energy, and a doubling of 
the scattering line width. If the Lorentzian line prohle is substituted in the analysis, the 
changes are slightly less signihcant. 

The observations of continuum cutoh energy relied almost entirely on the HEAO-1 
and RXTE measurements. The cutoh energy could not be measured with the many 
balloon experiments shown in Table 3 and Figure 3 because of atmospheric absorption 
below 25 keV. The MIR/HEXE experiment had an electronic threshold at 25 keV and was 
therefore insensitive to cutoh energy. The BeppoSAX ht was made with a more complicated 
continuum model with two inhection points near 20 keV, so could not be used. Einally, 
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errors were not stated for the Ginga observation. The most striking feature of the behavior 
of the cutoff energy is its apparent close proportionality to the centroid energy, shown in 
Figure 4. The dashed line displays proportionality with the best-£t factor of 2.03. The 
correlation coefficient for these data is 0.97, with a null (no correlation) probability of only 
1.1 • 10 -®. 

We conclude that a historical change has occurred in the values of at least two 
characteristics of the Her X-1 spectrum: the cyclotron scattering energy and the spectral 
cutoff energy of the power-law continuum. 


4. Discussion 

Models for the overall X-ray spectra and pulse shape from X-ray pulsars (Kirk et 
ah 1986; Meszaros & Nagel 1985) have had only moderate success in predicting the 
observed continuum spectral shapes. We concentrate therefore on the interpretation of the 
cyclotron line, especially in its bearing on the pulsar magnetic held. We discuss also the 
correlation of the continuum cutoff with centroid. The absence of sizeable change in the 
other spectral parameters should perhaps not be surprising given the likelihood that an 
important controlling factor in the spectral shape is the rate of energy release, and that this 
varied only over a factor of two in the joint RXTE and HEAO-1 data sets, for which we can 
speak most authoritatively. The stability of the parameters of the fluorescent iron emission 
at 6.5 keV and their similarity to the values obtained with OSO-8 (Becker et al. 1977) and 
BeppoSAX (Dal Fiume et al. 1998) indicate that the geometry of the reprocessing site, 
usually thought to be at the magnetopause (e.g. McCray et al. 1982), has not changed 
dramatically. 



16 


4.1. Cyclotron Line Centroid Variability with the RXTE Dataset 

The absorption feature is widely regarded to result from resonance scattering of 
electrons in a magnetic field (Triimper et al. 1978) of intensity ~ 3 x 10^^ Gauss (Voges 
et al. 1982), given the earlier measurements of 35 keV for the centroid. These hve RXTE 
observations give a new average value of 3.6 x 10^^ Gauss and can be used to set limits to 
the range of magnetic field strengths sampled. These RXTE data are consistent with no 
variability of line centroid to within 3% RMS. The held strength at the scattering site is 
therefore also stable to a precision of 3%. The magnetic moment can be determined from 
the inferred held strength, to be 1.8 x 10^'^ Gauss-cm^, assuming magnetic dipole structure, 
scattering at the polar surface of the neutron star and a 10 km radius. For scattering 
somewhere above the surface, this limit constrains the range for the average height of 
scattering to within 180 meters (90% conhdence) in a dipole held, and requires a somewhat 
higher magnetic moment. 

Next, we test the RXTE data against a more specihc process. There seems to be no 
general agreement on the structure of the radiating region of the polar cap of an X-ray 
pulsar. If a shock forms in the how above the pole BAK (1991) predict an accretion mound 
of sizeable extent. Other theorists, however, most recently Bulik et al. (1995), place the 
site of deceleration of the accretion how at the base of the neutron star atmosphere, with a 
scale height of less than a meter. The present results are consistent with the latter scenario, 
but they do not constrain the accretion mound scenario unless this is hundreds of meters 
or more in height. However, we note that if the observed line width of 5 keV or 12% is due 
largely to a mix of magnetic held strengths in a dipole held, then the scattering occurs over 
a range of the order of 1 km in height. 

BAK estimate an accretion mound height of 650 meters for Her X-1, with a direct 
proportionality of this height to mass transfer rate (luminosity). Although BAK state that 
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the radiation escapes from the mound at its base, Mihara et ah (1999) have had some 
success modeling centroid variability in several pulsars with intensity under the assumption 
that the scattering occurs at the top of the mound. Mihara et ah (1999) do not hud such 
a clear pattern, however, for Her X-1. As shown in Figure 5, the present RXTE data also 
rule out the BAK/Mihara centroid-intensity variability, to better than 3a conhdence. The 
RXTE observations span a 2-30 keV flux range of (4.5-8)xl0“® erg-cm“^-s“^, for which 
BAK/Mihara predict a change of 7 keV. In addition, there is no evidence in any of the 
RXTE observations for photoelectric absorption, therefore an actual correlation could not 
have been disguised by unrelated obscuration of the X-ray emitting region. 

The fairly tight RXTE constraints on variability of line centroid, interpreted as a limit 
on the range for height of scattering in a dipole magnetic field, give some preference for 
emission at the stellar surface. We note, however, that on the pulsation time scale the 
change of line centroid with pulse phase reported by Soong et ah (1990) could be produced 
in the BAK/Mihara picture if the average height of scattering changes with viewing angle, 
i.e., pulse phase. It has not yet been decided observationally, however, whether it is the 
centroid that changes more with pulse phase or the line prohle more generally. 


4.2. Variability on long Time Scales 

This new observation of a long-term change of cutoff energy closely proportional 
to the change of scattering energy suggests a common factor in producing both spectral 
features. This controlling factor, undoubted in the case of the scattering feature, is almost 
certainly the magnetic field strength. In a study of twelve X-ray pulsars with cyclotron 
lines Makishima et al. (1999) demonstrated a convincing correlation between the cutoff 
energy of the continuum and the line centroid. Their functional relation is a power law 
with index 1.4, thus cutoff ~ centroid^ "^. Could this population characteristic also apply 
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to changes in time with a single pulsar, as in this case? We hnd (Figure 4) that a power 
law relation with index 1.4 hts the RXTE and HEAO-1 data much worse than simple 
proportionality. However, we note that the data in Figure 4a of Makishima et al. (1999) 
seem reasonably consistent also with direct proportionality, so that there may indeed be 
some correspondence. 

These long-term spectral changes could result from a change of the height of scattering, 
although we can note that in this case the BAK/Mihara scenario would require HEAO-1 
fluxes two to three times what was observed, and thus can be ruled out. 

Brown & Bildsten (1998) have reconsidered the structure of the settling mound 
addressed earlier by BAK and have developed a model with a more compact structure, of 
the order of 100 meters in height. In calculating the hydrostatic pressure of the accreted 
settling matter in this structure, they obtain a pressure comparable to the magnetic 
pressure, which distorts and spreads the magnetic held lines, thus lowering the magnetic 
held intensity. They give no time scale for signihcant change of this structure. The observed 
limit of 5 years to the period of change may be applicable to their model. 

Ruderman et al. (1998) have a model of neutron stars in which motions of superhuid 
neutron vortices may alter the magnetic held in the core, which can result in movements of 
the stellar crust which in turn alter the surface magnetic held. The observed change of line 
centroid could thus in principle be related to a change of magnetic held in response to a 
change of structure. A structural change, in turn, would very likely produce a timing event. 
However, it is not known whether such a crustal change could produce a timing glitch large 
enough to be seen in the presence of accretion timing noise. 



5. Conclusions 


We have shown convincing evidence for a 23% historical change in the cyclotron line 
energy in the spectrum of Her X-1 through comparison of RXTE observations with earlier 
observations. For this purpose archival HEAO-1 spectra were re-analyzed, using modern 
tools and spectral forms for consistency. On the other hand, the hve RXTE observations, 
which span about two years, show no change, with a limit of about three percent. We also 
demonstrate a historical change of continuum cutoff energy which is closely proportional to 
the cyclotron centroid energy. This is the first good evidence to connect the cutoff energy, 
which is common to accreting pulsars, to the magnetic held. These RXTE data do not 
support a particular accretion mound scenario. 
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Appendix: Calibration of the HEAO-1 A4 LED detectors 

The two Low Energy Detectors (LED’s) of the High-Energy X-ray and Low-Energy 
Gamma-Ray Experiment (A4) aboard HEAO-1 were thin Nal scintillators of area 103 
cm^ each, energy range 12-180 keV, and energy resolution 25% FWHM at 60 keV. Many 
spectral results were published from these detectors, but only a very brief description of 
the energy calibration and monitoring procedures has been reported (Jung 1986). For 
energy calibration the LED detectors were provided with calibration sources whose 

events were recognized through a special pulse shape. But this calibration required special 
operating modes that were employed infrequently. Moreover, unexpected changes in the 
crystal optical interfaces following launch left the absolute energy calibration with these 
sources uncertain. Therefore an alternative procedure was devised which relied on the 
correct identihcation as well as monitoring of spectral lines in the instrument internal 
background. This background resulted largely from radioactive daughter isotopes produced 
by spallation interactions of cosmic rays and geomagnetically-trapped protons with the 
iodine nuclei in the detector. Most daughters are unstable and beta decay, producing a 
continuum in energy within the detector. But an alternative channel for positron decay 
is capture by the nucleus of an atomic electron in the K shell. This K-capture produces 
characteristic radiation whose total energy results from the fluorescent radiation of the 
electron cloud as it reconhgures as well as the nuclear deexcitation gamma ray, if any. 

A number of such lines were identihed in the LED detectors (Gruber et ah 1988). The 
most useful line for calibration and monitoring was from (half life 60 days), for which 
there was a gamma with energy 35 keV and a total of 32 keV from the electronic cascade. 
Because all the charge was collected by the electronics in less than the several microsecond 
pulse shaping time, the summed pulse corresponding to 67 keV energy loss was observed. 
The identihcation of this isotope and line in the A4 data was quite secure, because the 
60-day buildup of the activity at the beginning of the mission was clearly observed. 
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The most important complication for energy calibration for the A4 detectors was the 
energy-dependent efficiency of light production by the inorganic phosphor (Adams & Dams 
1970). Light production in the crystal lattice is stimulated by fast electrons resulting from 
the photon interaction. This process apparently has a maximum efficiency at electron 
energy of 10-20 keV (e.g. Zerby, Meyer & Murray 1961). This variation in light output 
with energy for a HEXTE detector is shown in Figure 4 of Wayne et al. (1998). The 
energy-dependent efficiency was carefully calibrated prior to launch of the HE A 0-1 using 
radioactive sources. This calibration curve was then used to calculate the light production 
for a decay, whose energy is shared among several electrons ranging from a few keV to 
tens of keV. The light production was calculated to be equivalent to light originating from 
a single photon of 62.7 keV, not 67 keV. Use of the former number successfully resolved 
a problem in the interpretation of LED spectra: a strong spurious 40 keV emission line 
which appeared in the spectra of the Crab, Cyg X-1, and the diffuse X-ray background! 
Reassignment of the calibration energy from 67 keV to 62.7 keV completely removed this 
artifact in all these spectra . From the sensitivity to this feature alone we estimate the 
absolute accuracy of the A4 LED calibration as 2%. Other aspects of the calibration were 
discussed most fully in several UCSD PhD dissertations: effective area, Nolan (1982); 
angular response, Knight (1981); gain variability, monitoring and correction, Jung (1986). 
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Fig. 1.— (Top panel) Fit with high-energy cutoff continuum and Gaussian absorptive line 
prohle to PGA and HEXTE data from 1996 July in units of detector counting rate versus 
energy. Except above 70 keV, where adjacent channels have been grouped for display, original 
detector energy channels are shown. These oversample the actual detector energy resolution 
by factors of three to eight. (Middle panel) Residuals to £t without absorption line, showing 
strong and correlated deviations, particularly at 40 keV, the line energy. (Bottom panel) 
Residuals, now acceptable, to fit which includes the absorption line. PGA errors include 
a 0.25% systematic error, which dominates the statistical errors. HEXTE markings, x for 
Gluster A and o for cluster B, have been omitted in the top panel for clarity. Error bars are 
one sigma. 

Fig. 2.— Spectral fit and residuals to HEAO-1 data of 1978 Feb 24. The three panels are as 
in Figure 1. The 1-16 keV data are from an Argon-hlled detector of the A2 experiment, and 
the 12-100 keV data marked with circles and crosses are from the two Low Energy Detectors 
of the A4 experiment. Errors (one sigma) are purely statistical. Original energy channels 
are shown, and these oversample the detector resolution by factors of two to four. Other 
than a negative feature at 4.5 keV, the fit is acceptable. The best-fit cyclotron energy is 33.2 
keV and the cutoff energy 15.9 keV. 

Fig. 3.— History of centroid energy of absorption line. Errors are one sigma, as estimated 
by the original investigators or in this paper for RXTE and HEAO-1. A pronounced change 
is evident following the Ginga measurement in 1990. Dashed lines indicate the averages for 
the earlier and later data. A similar and proportional change of the cutoff energy is also 
quite significant. The internal consistency of the RXTE data indicates that the errors have 
not been underestimated. Observations are listed by date and experiment in Table 3. 

Fig. 4.— The value of the centroid energy of the cyclotron scattering versus the cutoff energy 
of the spectral continuum, with 1-a errors. Gorrelation is formally highly significant, and 
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the dashed line shows the best-fit proportionality with constant 2.03. Errors are determined 
from the 68% joint conhdence regions. A power law dependence with index 1.4 (dot-dashed 
line), as was observed in a population study by Makishima et ah (1999), hts quite a bit 
worse. 

Fig. 5.— Measured line centroids versus 20-60 keV flux obtained with HEXTE-only datasets, 
with 1-a errors. The line shows expected dependence predicted by Mihara (1999) from the 
model of BAK (1991), in which the mean luminosity corresponds to a height of 1 km above 
the polar cap for the scattering. The predicted dependence is not conhrmed. 
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Table 1: Outline of Observations 


Observation 

Date 

Duration in sec 

Binary Phase Range! 

1 

1996 Feb 

5000 

.17-.23 

2 

1996 July 

28000 

.63-.93, .12-.46 

3 

1996 Oct 

5800 

.06-.08, .49-.52 

4 

1997 Sept 

20000 

.13-.38 

5 

1997 Nov 

18000 

.21-.40, .68-.74 


1 Zero phase is mid-eclipse, or superior conjunction 







Table 2: Results of Spectral Fitting, PCA <25 keV plus HEXTE Data 


date 

1996 Feb 

1996 July 

1996 Oct 

1997 Sept 

1997 Nov 

all 

2-30 keV flux, units 10 

keV cm” 

sec“^ 




flux 

7.84 

7.27 

5.33 

6.83 

4.44 

- 

power law with high-energy cut-off continuum 


indexl 

0.855(4)^ 

.862(4) 

.841(5) 

.901(4) 

0.852(4) 

.841(2) 

cutoff 

19.83(11) 

20.23(8) 

19.10(12) 

20.27(9) 

19.50(10) 

19.60(9) 

fold 

10.87(22) 

10.65(15) 

10.02(19) 

10.49(19) 

10.08(16) 

10.79(9) 

Fe line 







eqw 

.154(19) 

.190(18) 

.219(24) 

.140(16) 

.188(23) 

.171(9) 

centroid 

6.625(31) 

6.518(24) 

6.524(25) 

6.490(31) 

6.393(33) 

6.535(15) 

sigma 

. 31 ( 4 ) 

.37(3) 

.34(3) 

.30(4) 

.48(4) 

.36(2) 

scat, line 







depth 

0.80(5) 

0.74(3) 

0.63(3) 

0.66(5) 

.67(4) 

.75(2) 

centroid 

40.9(5) 

40.8(3) 

39.4(7) 

40.1(4) 

40.0(6) 

41.0(2) 

sigma 

4.5(5) 

4.9(3) 

3.8(7) 

4.7(4) 

3.7(6) 

4.9(2) 

Xr 

1.10 

0.99 

0.94 

0.98 

1.08 

4.56 


1 parameters have units keV except for index and depth, both dimensionless. 
^ 68% statistical errors in parentheses apply to the last signihcant £gure(s) 




























Table 3: Historical Observations of Parameters for Cut-off Continuum 


year 

Data Set 

r 

E, 

E/ 

T 

Er 

a 

+ 1900 



keV 

keV 


keV 

keV 

77.671 

BalloouHEXE 



12.2(9) 

3.2(6)t 

36.8(14) 

1.3 

78.151 

HEAO-1 


16.4(9) 

10.8(4) 

1.20(55) 

34.0(8) 

2.9(13) 

78.161 

HEAO-1 


16.6(10) 

11.1(5) 

1.48(98) 

33.9(8) 

2.5(14) 

78.164 

HEAO-1 


16.1(15) 

11.7(7) 

1.78(81) 

34.5(9) 

3.4(13) 

78.622 

HEAO-1 

.943(14) 

16.0(4) 

10.9(4) 

1.16(37) 

33.2(9) 

2.6(13) 

78.630 

HEAO-1 

.912(12) 

16.8(4) 

11.4(4) 

0.87(9) 

34.6(9) 

4.5(11) 

78.638 

HEAO-1 


16.9(5) 

9.8(4) 

0.68(17) 

33.2(9) 


80.726 

GRIS 



9.9(16) 


35.4(21) 

3.59(188) 

87.528 

HEXE 



10.6(10) 

3.0(18) 

33.8(18) 


87.618 

HEXE 



9.9(10) 

1.6(10) 

32.3(23) 


88.565 

HEXE 



10.8(20) 

4.5(34) 

35.0(21) 


90.070 

Ginga 

0.93 

19.2 

15.3 


36.8 

6.75 

94.370 

BATSE 





44(1) 

2.3(10) 

96.085 

RXTE 

.855(4) 

19.83(11) 

10.87(22) 

0.80(5) 

40.9(5) 

4.5(5) 

96.562 

BeppoSAX 




0.73(3) 

40.3(2) 

6.3(4) 

96.567 

RXTE 

.862(4) 

20.23(8) 

10.65(15) 

0.74(3) 

40.8(5) 

4.9(3) 

96.762 

RXTE 

.841(5) 

19.10(12) 

10.02(19) 

0.63(3) 

39.4(13) 

3.8(7) 

97.688 

RXTE 

.901(4) 

20.27(9) 

10.49(19) 

0.66(5) 

40.1(7) 

4.7(4) 

97.860 

RXTE 

.852(4) 

19.50(10) 

10.08(16) 

0.67(4) 

40.0(6) 

3.7(6) 


1 Equivalent width (keV) for subtractive Gaussian 
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Table 4: Average Values of Spectral Parameters in Table 3 



P 

E, 

E/ 

r 


a 



keV 

keV 


keV 

keV 

Pre-1991 

.927(8) 

16.52(17) 

10.84(20) 

0.90(11) 

33.67(25) 

3.32(34) 

Post-1991 

.864(10) 

19.89(21) 

10.40(16) 

0.71(3) 

41.02(40) 

4.83(36) 

change 

-.063(13) 

3.37(27) 

-0.44(25) 

-.19(11) 

7.36(47) 

1.51(49) 

ratio 

.932(13) 

1.204(18) 

.959(23) 

.789(102) 

1.218(15) 

1.455(184) 

(no change), dof 

210.7/7 

281.5/10 

35.2/15 

47.7/16 

831.6/20 

35.2/12 

(change), dof 

128.7/5 

97.1/8 

29.6/13 

41.5/14 

54.7/18 

28.1/10 

P null 

0.29 

0.014 

0.33 

0.38 

2.3 • 10“^^ 

.33 
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